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We compute the photon production of a QCD plasma at leading order in the 
strong coupling with a finite baryon chemical potential. Our approach starts from 
the real time formalism of finite temperature field theory. We identify the class of 
diagrams contributing at leading order when a finite chemical potential is added and 
resum them to perform a full treatment of the Landau-Pomeranchuk-Migdal (LPM) 
effect similar to the one performed by Arnold, Moore, and Yaffe at zero chemical 
potential. Our results show that the contribution of 2 i-> 3 and 3 4 2 processes 
grows as the chemical potential grows. 
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I. INTRODUCTION 



Heavy ion collision experiments at RHIC typically create a medium where the net baryon 
density is non- vanishing [HE]- As we enter a new era of precision measurements, it is 
therefore important to consider the effect of a non-vanishing baryon chemical potential on 
the thermal photon yield of the quark gluon plasma. Previously, the effect of non-zero 
chemical potential in 2 -h- 2 processes was studied in J3HH]. The complete leading order 
calculation of photon production including the effect of collinear enhancement in the 2 — > 3 
and 3—7-2 cases was first carried out by Arnold, Moore, and Yaffe (AMY) starting from 
first principles within the framework of quantum field theory at finite temperature [TJ. The 
analysis of [7J , however, is mostly carried out at zero chemical potential. The effect of adding 
a chemical potential was mentioned, but a detailed analysis was not fully performed. 

The chemical potential will modify the quark and gluon self-energies as well as change the 
statistical factors, thereby potentially modifying the power counting analysis of AMY. In 
this paper, we will precisely determine under which circumstances the power counting must 
be modified to account for the presence of a chemical potential. We also explore numerically 
the consequences of including a chemical potential in the photon production. 

As in the previous cases, the most convenient basis to work in when analyzing the para- 
metric sizes of diagrams is the Keldysh or r, a basis whereas the rates are most conveniently 
written in the usual 1, 2 basis. The addition of a finite chemical potential also changes the 
way in which we switch from the 1, 2 basis to the r, a basis, and we will carefully explain 
the required changes. 



The organization of this paper is as follows. In section II. 1 we briefly recall the structure 



of perturbation theory in the real time formalism and outline the basic formula giving photon 



production in terms of Feynman diagrams. In section II. 2 we show how to generalize a result 



of Heinz and Wang [8 J that allows one to express Green functions in the 1, 2 basis in terms 



of a reduced set of Green functions in the r, a basis. In section II. 3 we show how a finite 



chemical potential alters quark and gluon thermal masses. Then, in section III, we perform 



a power counting analysis to determine which diagrams contribute at leading order in the 



strong coupling g and resum these in section IV Finally, some numerical results are shown 
in section [V] 
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II. DIAGRAMMATIC APPROACH TO CALCULATING PHOTON 

PRODUCTION 

II.l. Perturbation Theory at Finite Temperature and Chemical Potential 

We start by briefly outlining the structure of real time perturbation theory at finite 
chemical potential. Since we are interested in a QCD plasma, our starting point is the QCD 
Lagrangian: 

c = E *M - M - 9^T a )iP f - \f%f? (ii.i) 
/ 

where the sum is over the Nf fermion flavors and the gauge group is SU(iV c ) with iV c = 30 
This Lagrangian has a conserved charge Q = Y2f J d 3 xi/j^(x)ipf(x) which is equal to the net 
fermion (quark) number. The density operator describing the grand-canonical ensemble is 
therefore e~ /3 ^~ At( ^- ) . In the imaginary time formalism, one can show that this changes the 
Matsubara frequencies from iu n to iu n + \x [91 [TUj. The chemical potential here is hence the 
quark chemical potential which is 1/3 of the baryon chemical potential. 

To connect the imaginary and real time formalisms, one defines the retarded and advanced 
propagators: 

G Ret (x) =ie(t)({iP(x)M0)}) Pttl (II.2) 
G Ad v(x) = -«0H)({V(aO,$(O)}W (II-3) 

The anticommutators above are replaced with commutators for gauge fields. By using the 
spectral representation of imaginary time propagators (for instance, see [91 [10]) one can 
show that the retarded propagator is obtained by analytically continuing the Matsubara 
propagator through the prescription iu n + /i i-> p° + it and the advanced propagator from 
the continuation iu n + fi (-> p° — ie, where p° is an arbitrary real number. 

In the real time formalism, one is forced to double the degrees of freedom because when 
constructing a generating functional from the partition function, the time integration contour 
has to traverse the real time axis once forward and backwards (Schwinger-Keldysh Closed 
Time Path) [11] . The first set of fields xpi, tpi, and corresponds to fields with a time 



1 Our metric convention is g^ v — diag(l, — 1, — 1, — 1). 
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argument on the forward directed part of the contour and conversely the set of fields ip2, 4>2, 
and corresponds to fields with time arguments on the backwards directed part. 

For computational purposes, it is sometimes more convenient to use another basis than 
the 1, 2 basis above. One such basis is the r, a or Keldysh basis [HI [12], defined as follows: 

ft /TT A \ 

ip r = (II.4) 

V?a = <fl ~ ¥2 (II-5) 

where p denotes any of our fields. By using the above algebraic relation and the spectral 
representation of propagators in the r, a basis, one can show that for fermions: 

G rr (P)=(±-n f ( J f > -ri)p(P) (II.6) 

G ar (P) = G Adv (P) (11.7) 
G ra (P) = G Rct (P) (11.8) 
G aa (P) = (II.9) 

The above allows us to compute real time propagators by analytically continuing imaginary 
time ones. From now on, we use capital letters for 4-momenta and lower case letters to 
denote the magnitude of the 3-momenta. For instance, P = (p°, p) in the above expressions 
and p — |p|. 

To tie this formalism to the problem of photon production, we recall the standard formula 
that relates the emissivity of the plasma to a Wightman current-current correlator [HI [TO] 

W = (M £ W(^"CO (n.io) 

v ' 1 1 a=l,2 

where we have defined the Wightman current-current correlator by: 

W^(K) = [ d'x e iK -*(f(V)f(x))^ (11.11) 



Wightman functions are most naturally expressed in the 1,2 formalism of the real time 
theory: 

, r _a EM d 3 k f d^P, f d*P 2 ((Pi\\+k) 2 +p\^ 

dT "-~^MJ wi (2,y 2 Pl{l ( Plll + k) ^- p ^ x 

x G 1122 (-P 1 , K + P l7 -K - P 2 , P 2 ) (11.12) 
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FIG. 1. The momentum assignments for computing the Wightman function. Arrows on fermion 
lines going into the shaded box correspond to insertions of the particle operator and those coming 
out of it corrrespond to insertions of the antiparticle operator. All momenta flow into the shaded 
box. 



where G1122 is a fermionic four-point function where two external vertices are of the '1' types 
and the others are of the '2' types. The "parallel" component of p is always defined relative 
to the fixed direction of the emitted photon momentum k. 

A schematic diagram for dY 1 is shown in figure [1} Our convention is that fermion momenta 
flow into the shaded box. Arrows going into the shaded box correspond to insertions of the 
particle operator ip, and arrows coming out of it correspond to antiparticle insertions -0. 



II. 2. Going from the 1,2 to the r,a Basis 



For power counting and actual computations, the r, a basis is more convenient than 
the 1, 2 basis. By using the algebraic relation between fields in the 1, 2 and r, a bases 
(c.f. Eq.(II.4[)), we can express £71122 as a linear combination of the sixteen possible four- 
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point functions in the r, a basis. The result is: 

1 1 1 

Cll22 G rr rr ~^G rarr -\- ~G arrr ~^G aarr 

- -Grrar ~ ifrrra + ^G rraa - -G rara 

~^G arra ~^G raar ~^G arar -\- ~G araa 

~^G aara -\- ~G raaa ~G aaar -\- G aaaa (11.13) 
o o o lb 

One can show that in general, G aa ___ a (l, 2, n) = 0, and hence we can immediately elim- 
inate one of the sixteen terms in the above - See for example [12J. In Ref. |Ej, Wang and 
Heinz showed that the remaining 15 terms can be re-expressed in terms of seven four-point 
functions and their complex conjugates. 

The authors of Ref. [8] only considered the case of a real scalar field. However, their 
derivation works equally well for fermions - See appendix A. Hence the following results 
from the reference still holds 

C1122 = OL\G aarT + OL2G aaar + a>3G aara + OL^G araa + OL^G raaa + OL^Garra + OC-jG arar 

+ 0iG* aarr + f3 2 G* aaar + f3 3 G* aara + (3<iG* araa + f3 5 G* raaa + (3 e G* arra + (3 7 G* arar (11.14) 

with the coefficients a« and composed of Fermi-Dirac distribution functions instead of 
Bose-Einstein distribution functions and G denoting the charge conjugate (ip O ip) of G. 
Here all G...(Q\, Q2, Q3, Qa) are functions of 4 momenta. Our convention is that the momenta 
Q2 and Q4 correspond to the insertion of ip and Q\ and Q3 correspond to the insertion of ?/\ 
This expression greatly reduces the number of diagrams we need to estimate. As a matter 
of fact, what emerges from our power counting analysis is that only G aarr contributes to 
photon production at leading order. This is because only this labeling gives rise to pinching 
poles as we will explain shortly. Thus, we only need to know the coefficients ot\ and (3\. 
These are given by 

ai = nf(qi+fi)nf(j^-fi) (11.15) 

This is proven in Appendix A. One can intuitively understand the signs of // in the above 
by recalling that nf(E — //) is the distribution function for particles while nf(E + ji) is the 
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distribution function for antiparticles. Since Qi and Q3 correspond to antiparticle insertions, 
they must be associated with the distribution rif(q® + /x), and conversely for Q 2 and Q 4 . 



As we will discuss in section III, all gluon exchange momenta must be soft at leading order 
in the strong coupling g. Therefore, we have that Pi ~ P 2 in figure [TJ Hence, Qx ~ —Qi 
and Q 2 — —Q3 and consequently, 0x ~ n/(g° + fi)nf(q 2 — ji) = ax- Further, one can verify 
that G aarr (-Px, K + Px, -K - P 2 , P 2 ) = G aarr (-Px, K + P u -K - P 2 , P 2 ) -See appendix B. 
Therefore, only the real part of G aarr is relevant to photon production at leading order. 



II. 3. Self-Energies at Finite Chemical Potentials 

When we compute the Wightman correlator to derive the photon production, we need 
to do a power counting analysis to identify all leading order diagrams. In this analysis, the 
appearance of "pinching poles" makes it crucial to resum the thermal self-energies into the 
quark and gluon propagators. Therefore, we need to know how the presence of a chemical 
potential affects the self-energies of quarks and gluons. This is well-known and dates back 
to the original paper by Braaten and Pisarski on hard thermal loops |13j . 

In order to keep our work self-contained, we review the derivation of self-energies with 
a full inclusion of a chemical potential in appendix C. In this section, we simply quote the 
final results. 

The gluon polarization tensor at finite chemical potential fi takes the form: 

IWP) = ml (-« + / (H.17) 

where in the above, m 2 D = g 2 ( j^f + ^jf~ + ^2 1 is the Debye mass, and = (1, v) with 
v = k/|k|. 

Similarly, the quark self-energy at finite chemical potential is: 

E T (P) = ^> (f + ^ / (11.18) 

The above two results show that the only effect of the chemical potential is to shift the 
dependence of self-energies on T 2 by a fi 2 term. As we will see, the consequence of this is 
that in carrying out our power counting analysis, we need not alter the arguments of Arnold, 
Moore, and Yaffe as long as /i < 0{T). 
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FIG. 2. The pair of fermionic propagators carry momenta which differ by K. We assume that K 
is nearly collinear with P and that the other parts of the diagram are essentially constant over the 
region of width g 2 T centered around the pinching poles of the integrand. 

III. POWER COUNTING WITH A FINITE CHEMICAL POTENTIAL 

When evaluating the Wightman function, there are two regions of the spatial part of the 
loop momentum integration that are of interest. 

• The non-collinear region: Pi,± and p2,± are both 0(T). 

• The near-collinear region: pi ^ or p2,j_ is O(gT) or less. 

The non-collinear region corresponds to the contributions of the basic 2 i— > 2 processes 
treated by Kapusta et al. [H] and Baier et al. [15] . The near-collinear region will have 
both pinching pole and near-collinear enhancements. It corresponds to the contribution 
of the bremsstrahlung and inelastic pair annihilation processes modified by the Landau- 
Pomeranchuk-Migdal (LPM) effect. 

We briefly review how pinching pole enhancements arise. As will be shortly shown, the 
leading order diagrams for photon radiation all contain a pair of fermionic propagators as 
shown in figure [2j The spatial part p of the loop momentum P is assumed to be nearly 
collinear with k. In other words, p^ ~ gT. Given that we impose the collinearity of p and 
k, we may focus just on the frequency integral 
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FIG. 3. The locations of the four poles of the frequency integrand when p and k are nearly aligned. 
The poles above the real line belong to G a dv(K + P) and the ones below belong to G TCt (P). The 
integration contour is along the real axis. We see that two of the poles nearly coincide at pit . These 
poles are said to "pinch" the integration contour. 



where we have absorbed the Dirac matrix structures into the vertices. In Eq.(III.l), T p = 
Im(E(P))/(2£'p) is the decay width of the quark generated by the imaginary part of the self 
energy. To leading order in g, it can be replaced by its asymptotic value T = lim p _ s>00 (r p ) 

(Bsnzi). 

The poles of the integrand are at the following locations: p° = ±E P — |r and p° = 
—k° ± -Ep+k + |r. It is straightforward to see that -E p +k ~ P\\ + k when p, k are 0(T) 
and collinear. Hence the two pole positions p° = E p — iT/2 and p° = — k° + E v+ ^ + iT /2 
almost coincide at pu = p ■ k although they are located on opposite halves of the contour 
(c.f. figure [3]). 

We use the residue theorem to close the contour and pick up the contributions of the 
poles located below the real axis. For the case where p\\ > 0, this gives: 



dp 
2^ 



G Adv (K + P)G Rct (P) 



2E p (k + E P - E p+k - iT)(k + E p + E p+k - iT) 



(III.2) 
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FIG. 4. All ladder diagrams contribute at leading order when the gluon exchange momenta are 
soft. 

When pii < 0, the factors in the denominators should be replaced by the general expression 

E p - E p+k + k -» E p signal, ) - E p+k sign^, + k) + k = 5E (III.3) 
E p + E p+k + k - iT -> E p sign(pn) + E p+k sign(p|| + fc) + fc ~ 2(p|| + fc) (III.4) 

Therefore, the frequency integral is approximately: 

dn° 1 

> (K+p)WP) = aTM (IIL5) 

But T and are both of order g 2 T, while p and k are hard. Hence, we get a l/g 2 T 3 en- 
hancement from the frequency integral. These enhancements make a large class of diagrams 
contribute to leading order even though we would naively expect them to be subleading. 

Suppose that in the above demonstration we had not taken the integrand to be the 
product of an advanced and a retarded propagator but rather J GAdv(P)GAdv(P + K) or 
/ G^ ct (P)G^ ct (P + K). Then all poles would have been on one side of the frequency in- 
tegration contour only and we could have closed it on the side with no poles. Since the 
contribution from great circles at infinity vanishes, this means that the integral would van- 
ish. Hence pinching pole enhancements only arise when we have a retarded and an advanced 
propagator. This is the reason why only G aarr contributes to the leading order photon pro- 



duction. Due to the fact that there is no aa propagator (c.f. Eq.(II.9)) and that an interacting 
vertex must contain an odd number of a-fields [12], any other labelling leads to reduction of 
the number of pinching poles. 

We will illustrate general power counting arguments by estimating the size of the ladder 
diagram shown in figure |4j We start with the ladder diagram that has only a single rung 
(figure [5]). Notice that the assignment of the r, a indices ensures we have a retarded and an 
advanced propagator in each pair of quark propagators. To get a pinching pole enhancement, 
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FIG. 5. A ladder diagram with one rung only. It has two pinching pole enhancements and a further 
enhancement coming from the size of the soft gluon rr propagator. 

we need that the quark with momentum P be nearly collinear with the emitted photon. 
Further, if the gluon carries a soft exchange momentum Q ~ gT, then it cannot disturb the 
collinearity of the quark with the emitted photon. Therefore, the second pair of propagators 
also gives a pinching pole enhancement. So we have a g~ 2 x g~ 2 enhancement. A soft gluon 
rr propagator is of order g~ 3 T~ 2 because of the Bose-Einstein factor, provided \i < 0(T). 
However, a soft gluon also brings in a phase space suppression of order g 2 which cancels this 
enhancement. 

Gathering all powers, we have for suppressions: 

• g 2 from the integral over P± which is 0(g 2 ) since this is the near collinear regime. 

• g 3 from the soft Q spatial integral over Q. 

• g 2 from two gluon exchange vertices. 

• g 2 from the P± ■ {P+Q)± factor arizing from the contraction with the external photons, 
and for enhancements: 

• g~ 2 x g~ 2 from two pinching pole enhancements. 

• g~ 3 from the soft gluon propagator. 

Adding up all powers of g, we get that this two-loop diagram is of order g 2 , the same as the 
non-collinear one-loop diagram. 

The above analysis easily extends to ladder diagrams with an arbitrary number of rungs. 
Each additional rung brings with it one more pinching pole enhancement, giving a g~ 2 , and 
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FIG. 6. The resummed vertex T>^ sums all ladder diagrams with one photon vertex removed. All 
external legs are amputated. 

a soft gluon propagator of order g~ 3 . It also brings in a g 3 suppression from the spatial 
integral over the new soft gluon momentum and a g 2 suppression from the two additional 
gauge boson exchange vertices. Adding these up, we get a net contribution of g°, so that at 
least all ladder diagrams contribute to leading order. 

One ought to make an important remark at this point. Since the thermal mass of the 
gluon goes roughly as T 2 + fi 2 /2ir 2 , the spectral function can be of higher order than g~ 2 if 
fi is of order T/g. In this case, the size of the soft rr gluon propagator fails to cancel the 
g 3 phase space suppression that we get from a soft gluon loop momentum. Consequently, 
all ladder diagrams are sub leading when /i = 0(T/g), and it becomes unnecessary to resum 
the contributions of all ladder diagrams. 

The rest of the power counting analysis aims at identifying which combinations of exter- 
nal r, a indices contribute at leading order, showing that all other diagram topologies are 
subleading, and also proving that gauge boson momenta of order higher or less than gT 
also give subleading contributions. As we have argued, as long as /i T/g, the size of the 
propagators comprising each diagram is parametrically the same as in the analysis of [7]. 
We will not pursue it here in more detail. 



IV. RESUMMATION OF LADDER DIAGRAMS 



In the regime where \x < 0(T), the power counting analysis of [7j is unchanged. There- 
fore, the conclusion that one needs to resum ladder diagrams (and only these) to get the 
photon emissivity to leading order is also valid. Hence, we can apply the same resummation 
procedure as [7J. Here we briefly outline that resummation procedure. 

Graphically, the sum we need to evaluate is illustrated in figure [6j We have defined the 
photon resummed vertex as the sum of all ladder diagrams with the leftmost photon 
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F{P-K)= M(P,Q,K) 



P + K 



P + K 

FIG. 7. The graphical operator T adds a pair of propagators evaluated at the location of the 
pinching poles. The rung Ai adds a soft gauge boson propagator. 

vertex and pair of quark propagators removed. From figure [6j one sees that each lad- 
der is constructed from the previous one by concatenating a gluon rung attached to two 
quark-gluon vertices and a pair of quark propagators to the previous one, where the quark 
propagators must be evaluated at the location of pinching poles. The ladder diagrams form 
a geometric series. 

= X M + MF1„ + MTMTT^ + MTMTM JT„ + ... (IV.l) 

The graphical operator T adds the pinching pole pair of propagators. In the above equation, 
each T contributes a pair of pinching poles: 



T{P; K) = G Adv (K + P)G Rct (P) 
-1 1 



4tt5 [2p° + k°-E p signal) - £ p+k sign^y + k)] (IV.2) 



4pn On + k) T + i5E 
The operator M. adds the rung as illustrated in figure [7] 

M(P,Q } K) = tg 2 C F 4 Pl] ( Pl] + k)K^k u G^{Q) (IV.3) 

Finally, X M is the bare photon vertex on the far end of the ladder diagram as illustrated in 
figure [8j 

The resummed vertex T>^ then satisfies the following integral equation: 



V tM (P,K)=l f ,(P,K)+ [ M{P,Q,K)F(P + Q,K)V^P + Q,K) 

Jq 



(IV.4) 



After defining the variable / M (p, k) = 4p||(p|| + k) J ^^(P, K)V^(P 1 K), the integral equa- 



tion (IV.4) becomes [7]: 

d 3 



2P» + K" = i6EfP(p, k) + J 7^C(q) [/" (p, k) - r(p + q ± , k)] (IV.5) 
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FIG. 8. Bare quark-photon vertex 1^. 



After applying a sum rule, the collision kernel C(q) is found to have the simple form [18] 

27r v_±/ q 2 ± q 2 ± +m 2 D 
with p and k strictly colinear. The term SE is the difference between the locations of the 



pinching poles of the quark propagators, as seen in Eq.(III.3). Through the energies of the 
quarks, it includes their asymptotic thermal masses m^. Therefore, the chemical potential 
enters through both and m^. Recall that 

2 2 fN f T 2 N C T 2 N f fi 2 \ , > 



4 V 7T 2 



mi = ^^ L \T l + ^\ (IV.8) 



Finally, the contribution of the 2 y 3 and 3 t— >■ 2 processes to the Wightman correlator 
(and thus to the photon emissivity) is given by: 

W% M {K) = J 0^ iP l^llf MP° ~ n/(p° ~ ^)]x 

x Z"(P, AT)Re[.F(P, K)V U (P, K)\ (IV.9) 



V. NUMERICAL CALCULATIONS 



Following a method developed by Aurenche et al. [19], we can solve equation (IV.4) 



numerically by solving the differential equation in impact parameter space. As in Ref . [20J , 
we decompose the photon emission rate as follows 

1 



\n(T/m c 



ln(2fc/T) + C 2 » 2 (k/T) + C hrem+ann (k/T) 



(V.l) 



where 



A(k) = 2a E M dp ^ Q: 



2 \ m oc 



n f (k) 



(V.2) 



r=u,o 
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FIG. 9. Combined q and q bremsstrahlung contribution to the photon emission. The ratio against 
the // = case increases with increasing /J./T. 

Here, k is the magnitude of the emitted photon's momentum and n/ is the Fermi-Dirac 
factor. 

Since the 2^2 part of the spectrum has been already calculated [3H6|, we only plot 



Vb+ a (k) = A(k)C\> Tcm+ann (k /T) in figure 12 The plots show the contribution of the 2^3 



and 3 (—7- 2 processes with the full treatment of the LPM effect. The temperature of the 



plasma was taken to be 250MeV. We see from figure 12 that as the ratio pi/T increases, 
the photon production rate also increases. This finding is consistent with the behavior of 
the 2 -f->- 2 emission rates. 

One of the reasons for this enhancement turned out to be just the statistical factors in 



Eq.(IV.9). These factors represent 3 different processes depending on the sign of p° and the 
relative sizes of \p°\ and k° = k > 0. To illustrate the effect of non-zero \x in each of the 



three processes, consider the following ratio of the statistical factors in Eq.(IV.9) and the 

same statistical factors with \x = 0. 

= n f ip + k -fi)(l-n f {p -fi)) 

n f (p + k )(l-n f (p )) 1 J 

When p° > 0, the underlying physical process is the bremsstrahlung from the quarks. In this 
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FIG. 10. Pair annihilation contribution. The ratio against the = case decreases with increasing 
M/T. 

case, the ratio is mostly greater than 1. Hence the rate is enhanced. This reflects the fact that 
a positive chemical potential enhances the number of quarks more than the Pauli-blocking 
factor reduces the emission rate. When p° < 0, (1 — rif(p° — //)) = rif(\p°\ + //) becomes the 
anti-quark phase space density. For p° < and k < \p°\, the factor 7i/(— \p°\ + k — fi) = 
1 — rif ( \p° | — k + /i) represents the Pauli-blocking factor for the anti-quark bremsstrahlung. In 
this case, a positive chemical potential reduces the phase density of anti-quarks but enhances 
the Pauli-blocking factor. But since Pauli-blocking can never be enhanced above 1, the effect 
is to reduce the photon emission rate. The combined effect of quark bremsstrahlung and 
anti-quark bremsstrahlung is a net enhancement as shown in figure |9j 

Forp < and k > \p°\, the factor n/(— \p°\ +k — jj) represents the density of quarks that 
can annihilate with anti-quarks to produce a photon with energy k. In this case r is always 
smaller than 1, reflecting the fact that the annihilation process is necessarily controlled by 



the lesser number of anti-quarks. In figure 10 , we show the contribution of the annihilation 



process to the photon emission rate for various values of 

Among these three effects, the enhancement of the quark-bremsstrahlung dominates be- 
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FIG. 11. For [i of the order of the temperature T of the plasma, the photon emission rate increases 
by about 10% relative to the u = value. 



cause it increases much faster than the reductions in the annihilation contribution and the 
anti-quark bremsstrahlung. Therefore, overall, the effect of having fi > is to enhance the 



photon production at the same temperature as shown in figures 11 and 12 This trend was 
also observed in studies of 2 -h- 2 processes jl] . In one previous study of 2 -h- 2 processes [3] , 
it was found that increasing the chemical potential decreases photon production. However, 
this was for constant energy density instead of constant temperature. 

Phenomenologically, the baryon chemical potential at RHIC is about 30MeV and at 
SPS it is about 240 MeV[2T]. The quark chemical potential is therefore about lOMeV and 
80MeV for RHIC and SPS, respectively. With T = O(200 Me V), the enhancement at RHIC 
is negligibly small and it will be modest at SPS, no more than 5 %. 



VI. CONCLUSION 



In this paper, we have explored the effect of non-zero baryon chemical potential in thermal 
photon production. Following the zero chemical potential study [7j, we have computed pho- 
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FIG. 12. Plots of the ratios Vb+ a {^) / v b+a{^ = 0) for < /z < 10T. As (j,/T increases, so does the 
ratio. 



ton production from non-Abelian plasmas at leading order by resumming ladder diagrams 
to fully incorporate the LPM effect. After a careful analysis, we have found that as long 
as fi < 0(T), the formulation in [7] is still valid with appropriate changes in the statistical 
factors, thermal masses and the Debye mass. However, when \x = 0(T/g), resummation 
of ladder diagrams is no longer necessary because thermal masses become 0(T) instead of 
O(gT). Hence, inverse powers of thermal masses no longer provide enhancements. 

Numerically, it is found that the inclusion of a chemical potential up to fi ~ T enhances 
the photon emission rate moderately, up to 30 % when /i = T. This trend is also valid for 
hard photons from the 2 <H- 2 process [1]. Since the quark chemical potential is relatively 
small compared to typical QGP temperatures at RHIC and the LHC, we expect a relatively 
small effect from a finite /x on thermal photon production at RHIC and the LHC although 
it could be significant at SPS energies and also for the lower energy runs at RHIC. 
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Appendix A: KMS Condition with Finite fi 



In this appendix, our goal is to prove equation (11.14). This equation was proven in 
reference [5] for the case of scalar fields. Their argument exploits the KMS condition to 
derive a system of equations that allows them to solve for G rrrr , G rrra , G rrar , G rarr , G arrr , 
G r aari G raraj and G rraa in terms of G aarrj G aaarj G aaraj G araaj G raaaj G arraj G arar and their 
complex conjugates. Their procedure also goes through for Dirac fermions in the grand 
canonical ensemble provided we modify the usual KMS condition to include a chemical 
potential. 

Consider first the commutation relation between Dirac field operators and the conserved 
charge Q: 

[Q,ifj\ = --0 e^i[> = e-^ePri 
[Q, $ = ^ e?i4$ = e ^e^ (A. 1) 



For Green functions in the 1,2 basis, these commutation relations allow us to relate Green 
functions whose 1, 2 index assignments are opposite each other. Indeed, consider: 

G 2 ..2i..i(x 1 , ...,x n ) = (f(ip^(xi)...i/}^(x r )t/}(x r+ i)...i/}(x m ))x 

x T(ij\x m+l )..4\x l )ij(x l+l )..4(x n )))^ (A.2) 

Note that the 2 index refers to the space-time points X\, x m and the 1 index refers to the 
points x m+ i, ...,x n . We have also omitted spinor indices in the above for simplicity. Using 
the fact that e@ H is a time evolution operator in imaginary time (ip(t + i(3) = e~^ H ip(t)e^ H ), 
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we obtain: 

G2..21..1OC1, ■••) x n) — 

(f{e p6 ^{x\ + i/3)...^{x° r + t(3)^(x° r+1 + i{3)...i){x° m + 0)e~ p6 )x 
x T{^{x m+1 )...^{x l )^(x l+1 )...^{x n )))^ (A.3) 

The operators e ±l3H can be pulled outside of the reversed time-ordering symbol. Therefore, 
when writing out the thermal average explicitly in terms of the density operator e~P( H ~^ 
, we obtain: 



G2..21..1O1, ••• ; 3*n ) 

1 

z 

x e-^T(^(x m+1 )...^(x l )^(x l+1 )...ij(x n ))) (A.4) 



We want to commute the operators inside the reverse time-ordering operator past e 13 ^ 
and use the cyclicity of the trace to take them to the right of the operators inside the time- 



ordering symbol. Equation (A.l) shows that in doing this, we pick up a factor of when 
we are commuting a i]r operator and a factor of e~^ M for a ip operator. Let us define the 
symbol Oi to be equal to 1 if the space-time point X{ has a ijy operator insertion at it, or 
— 1 if it has a ip insertion. With this notation, we obtain the following generalization of the 
KMS boundary condition to n-point functions: 

G*2..21..l(2 ; l, %n) — 

L e ^i=^Tr(e-^ A -^T^\x m+1 )^ 

x f^j\x\ + i/3)...^(x° r + i/3)iP(x r+1 + i/3)...iP(x° m + 0))) (A.5) 
Finally, we use the well-known formula (j){t + a) = e adt (j)(t) to obtain: 



G2..21..1 (^l > x n ) — 

e /3E l|a ^ 2 (^ + ^,) (T( ^ t(Xm+l) .^t (xi) ^ (xm) .^ (Xn))x 

x f (^(x?)...^^ )^^ ^)...^^)))^ (A.6) 

We are allowed to take the e^'i°* =2 operator outside of the time-ordering symbol 
provided we take time-ordering to be given by the T* prescription. That is, the time- 
ordered n-point function is defined to be the path integral of the product of the n operators 



21 



weighted by the exponential of i times the action. For a discussion of this prescription, see 
for example [22j. 

We define the "tilde conjugate" G of a Green function G to have same 1, 2 index assign- 
ment but with the time-ordering symbols reversed. Therefore, we can express our previous 
result as: 

G 2 .. 21 ..i(x l5 ...,x n ) = e^^i^ ( ^ + ^ ) G 1 .. 1 2.. 2 (x 1 ,...,x n ) (A.7) 



In [8], Wang and Heinz use the fact that in momentum space, tilde conjugation is equiv- 
alent to complex conjugation. As this is a key fact, we proceed to extend it to fermion fields 
at finite chemical potential. For notational simplicity, we will focus on the case where there 
are four fermion operators two of which bear "1" indices, and the other two of which bear 
"2" indices. Although this is not necessary for this discussion, we have also included the 



7 matrices which appear in the current- current correlator 11.11 Extensions to other cases 
should be obvious. 



i 



G 221 i(Qi,Q 2 ,Q3,Q4) = / Y[d 4 Xi e iQi - Xi {f $(x 1 )>f1>(x 2 ))T$(x 3 )>fil>M))f ) , ll (A.8) 



i=i 



Taking a complex conjugation and changing X{ — > —Xi, we obtain: 



GWQi,Q 2 ,Q3,Q4) = / n^ eiQ! ' Xl < T ^(- X4 )^(- x 3))T(^(-x 2 ) 7 ^(-x 1 )))^ 

(A.9) 

Next, we make use of CPT invar iance. Recall that the Dirac bilinear ip^^ip transforms as 
follows under the anti- unitary CPT transformation B: 

e^(x) 7 ^(l/)0 _1 = -^(-yh^i-x) (A.10) 



Then, using anti-unitarity, thermal expectations of a product of fields can be related as 
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follows: 



(r(^(-x 4 )7^(-^))^(-x 2 )7^(-^i))>^ 

^^(n|e^-^ ) T(^(-a;4)7^(-^3))T(^(-a;2)7^(-^i))h} 



n 



= \ ^(en|9e-^-^T(^(-a;4)7^(-^3))T(^(-a: 2 )7^(-^))e- 1 |0n) 

n 

= ^(n1e-^ + ^T(6^(-a:4)7>(-^)e- 1 )T(0^(-^)7^(-^)e- 1 )|n / ) 



n' 



= ^(n>-^ + ^)T(^(x 3 )7^(^))^(x 1 )7^(^))|n / ) 

n' 

= (T(^(x3)7^M)T(^(x 1 )7^(^))> / 9,^ (A.11) 

We have used CPT invariance in commuting 6 past e~^ H and relabeling the eigenstates 
\n'} = <d\n). Further, note that the sign of /x changes as we commute 6 past e 13 ^. Combining 
the above with our previous result, we obtain: 

r 4 

G* 2211 (Q 1 ,Q2,Qs,Q4) = / n^ eiQ< " < W(-^)7^(-^3))T(i/;(-x 2 )7^(-^i)))^ 



i=i 

4 



i=i 

= G 2211 (Q 1: Q 2 ,Q 3 ,Q 4 ) (A.12) 

where (5 denotes the charge conjugate ^) of C Recapitulating, we have the two 

equations: 

G 2 ..2i..i(K 1 , #„) = e^K-^^..^..^, K n ) 
G 2 .. 21 .. 1 (K 1 ,...,K n )=Gl, 21 , A (K 1 ,...,K n ) (A.13) 

The above can be taken as the starting point to reproduce the derivation of Heinz and Wang. 
However, the latter also uses the energy conservation condition £\ k® = and the identity 
n b{ki) + n b{—ki) = —1. In order for the derivation of Heinz and Wang to proceed in the 
same manner, we need suitable generalizations of these to treat fermions at finite chemical 
potential. By charge conservation, we must have that Y17=i 0i — - In other words, the 
Green function in question must have an equal number of ift and ip insertions to not vanish. 
Therefore, we can define £f = + <Tj/i and maintain "energy conservation" Y^=i M = 0- 
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Also, by performing the customary replacement fib(k^) i— > — n/(Af), we preserve the relation 
n f>(X°) + n b(— k®) = —1 since nf(E) + nf{—E) = 1 in general. 

The only relations among the k^s that Heinz and Wang use is energy conservation, and 
the only property of the Bose-Einstein distributions they use is n^k®) + n&(— k®) = — 1. 
This is because after they have used the KMS condition like we have done above, their work 
consists of algebraic manipulations such as solving large systems of equations involving the 
Bose-Einstein distribution. In their work, it is possible to treat G and G* as independent 
variables rather than considering the real and imaginary parts of G. Hence, it causes no 



harm to their derivation to replace G* by G* as is required by (A. 13). We conclude that we 
could in principle perform the exact same manipulations as Heinz and Wang by relabeling 
every n^kf) as —nf(k^). The net result is that we make the same replacements in their final 
results, bearing in mind that ip^ insertions come with a — Uf(k° + y) and if) insertions come 
with —nf(k° — fi). 

In reference jS] , Heinz and Wang obtain 

«i = n h {ql)n b {q° 2 ) (A.14) 



and 



Pi = -(1 + n 6 (g 3 ))(l + n M \ + Mq o )+nbiq o ) 



(A.15) 



By making the replacements we have prescribed above, we get: 

"i = n f((ii + v) n f(Q2 -y) ( A - 16 ) 



and 



as required. 
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Appendix B: A Simple Identity Involving Charge Conjugates 

In this appendix, we want to show that G aarr (Ki, K 2 , K 3 , K^) = G aarr (Ki, K 2 , K 3 , K 4 ) 
It is sufficient to prove the corresponding identity in position space. 

G aarr {x 1 ,X 2 ,X 3 ,X 4 ) = (Tc(^aMlp a (x 2 )^rM'^r(^)))fi^ 
G aarr (x 1 ,X 2 ,X3,X4) = (T C {'ip a {x2)i> a {xi)i) r (x4)i) r (x 3 )))^- li 

= ^( n l^ (M)T cW^)i(ll)^Wllv(l3))|n> 
n 

= \ ^(Cn|C , e-«^)T c (^ a (x 2 )^(x 1 )^ r (^ 4 )^r(^3))C' t |C'n) 

n 

= \ ^(^|e-' 3 ^-^CT c (^(x 2 )^ a (x 1 )^(x 4 )^(a; 3 ))C t |n / ) 



^^(n>-^ # -^)T c (a4( a ; 2 )V' a (^i)C'tC'^( a ; 4 )^(a;3)C't)|n / ) 



n' 



Z 

n' 

— G aarr (xi, x 2 , x 3 , X4) 

Appendix C: Computation of Self-Energies at Finite Chemical Potentials 

In this appendix, we closely follow the zero \x analysis in [23J to compute self-energies in 
the hard thermal loop approximation with a finite chemical potential. 

Consider first the gluon self-energy 11^. The color structure of this tensor is trivial, so 
that 11^ = 5 ab Il^ u . We have four diagrams to consider: the quark loop, the gluon loop, 
the gluon tadpole, and the gluon ghost. Only the quark loop is affected by the presence of 
a chemical potential and hence we make it the focus of our attention. With all momenta 



labeled, this diagram is shown in figure 13 



Our strategy is to evaluate this diagram in the imaginary-time, invert it to get the gluon 
propagator, and then analytically continue the result to retarded frequencies to get the 
real-time ra propagator from which all other propagators in the r, a basis may be obtained. 
Proceeding, we first get: 

IMP*) = g 2 Tr(T a T a ) J {dK E } Tr{ ltM Sf\K E ) lv Sf\Q E )) (C.l) 
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In the above, the integration measure is {dK E } = T ^^l_ oc / (dk), where (dk) = is the 
spatial momentum integration measure. Since we have factored out the color structure of 
the gluon self-energy, there is no sum over a in the Tr(T a T a ) factor. Further, Ke = (iu r , k) 
is the fermionic loop momentum and Qe = Ke — Pe = (iw r ~ ^n, k — p) is the difference 
between the momentum of the gluon (P# = (iu n ,p)) and the virtual quark. The frequency 
iu n = 2imT is bosonic, and the frequency iu r = 27r(r+ |)T is fermionic. The color matrices 
T a are the generators of SU(3) in the fundamental representation. In this representation, 
the group factor Tr(T a T a ) is simply |. S F ^ denotes the fermion propagator time evolved 
using the "thermodynamics Hamiltonian" H — fiQ since this is what satisfies the periodic 
boundary condition (see Appendix A for KMS conditions at finite //). 

The cleanest way to compute the above is to use the spectral representation 

Sf\iu r ,k) = dre^ T (i- nf (k -fi))te-^ k °-^p (K) (C.2) 

Jo J-oo 271 

Indeed, this expansion of the propagator makes the frequency sum trivial. After performing 
the spin trace and employing standard identities involving the Fermi-Dirac distribution, we 
obtain: 

n^K.p) = 2?N, J (dk) J^-J ^ Po (K)p (Q)x 

x (K.Q V + Q,K U - g, v K ■ Q) W/(fc °~ & 7 W/(g ° - ^ (C.3) 
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For the spatial component of the self-energy tensor, the result is 

U i:j (iuj, p) = 2g 2 N f JWj^Jf {hq 3 + q&j + M^k^q - k • q)) 

'rif(E k + p)- n/(£ q + p) n f (E k - p) - n f (Eq - p) 



iu + Ek — Eq iu — -Ek + Eq 

+ (kiqj + qikj - SijiE^E^ + k • q)) 

1 - n f (E k + p) - n f (E CL - p) 1 - n f (E q + p) - n f (E k - p) 



iu + E k + 



iu — E k — E c 



(C.4) 



using po{K) = -^r (5(k° — k) — 5(k° + k)). The above expression holds for arbitrary external 
momenta. Henceforth, we only consider the hard thermal loop approximation so that p ^ k 
after the analytic continuation iu h-> p° + ie has been performed. It turns out that the 
self-energies in the hard thermal loop approximation are also valid for arbitrary external 
momenta P at the one loop order [161 Ej- Hence, we will freely use our results as the 
self-energies of the gluons and quarks when estimating the size of their propagators. 
Using the fact that p <C k, we obtain: 

v • p d 



U ij (p°,p) = -g 2 N f J(dk) 



ViV 



p° 



v • p + ie dk 



(n f (k - p) + n f (k + p)) 



+ (viVj - Sij)-(nf(k - p) + ra/ (A; + pi}) 
where v = k/|k|. After an integration by parts and an application of the result 

I l \ ^2rp2 ,,2 



(dk)k 



+ 



e 0(k+tM) _|_ I e P(k-n) _|_ I 



7T T /T 
h — 

6 2 



we get: 



U tJ (p°,p) = -g 2 N f 



2 



TT^T 2 

~~ 6~ 



dQ 



v^j - 5ij - 2(r,r/) 



v ■ p 



v • p + ie 



(C.5) 



(C.6) 



(C.7) 



By symmetry, we have J dflViVj = | J dQ5ij. Thus, we finally obtain the hard thermal loop 
part of the gluon polarization tensor. 



U t ,(P)=g 2 N 



2 



^T 2 
~ 6~ 



dQ 



p°ViVj 



An 3 p° — v ■ p + ie 



(C.8) 



By going through the same analysis, we can evaluate the other components of the gluon 
polarization tensor. If we include the contributions of the three other loops, we finally have: 
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Q 




p K 

FIG. 14. The chemical potential appears explicitly in the quark propagators, whence we need to 
resum the self-energy of the quark. 

where in the above, m 2 D = g 2 ^ — + + is the Debye mass, and = (1, v). 

The contribution of the quark loop to the gluon self-energy corresponds to the screening 
of the strong interaction by quarks in the medium. It is therefore natural to expect that the 
chemical potential of the quarks will have an influence on the repeated scattering events that 
occur during the photon emission process. The chemical potential also appears explicitly in 
the quark propagator and thus affects the thermal mass of the quarks. This will affect the 
integral equation that we will derive for photon production, so we turn to computing the 



self-energy of the quarks. The relevant diagram is shown in figure 14 

is: 



The expression corresponding to the diagram in figure 14 



E{P E ) = g 2 T a T a j {dQ E } llM S'P(P E - Q E )i v G^{Q E ) (C.IO) 

In the above, f {dQ E } = T Yl^L-oo I C^) denotes a sum over bosonic frequencies followed 
by an integration over the spatial part of Q. As usual, T a T a = C*2(-F)l = |l where 1 is the 
identity matrix in SU(3). In the Coulomb gauge, we find: 



S(iw n ,p) = -g 2 C 2 (F) J {dQ}(ry S , M{K E )'Y Q 



1 f dq° p L (g°,q) 



q 2 J 2tt q° — iu n 



where we have used the spectral representation of the longitudinal and transverse part of the 
Coulomb gauge gluon propagator following [23]. Here g« = qi/q. The above shows explicitly 
that we have three contributions to the quark self energy: E(P) = Sc(-P) + ^l(P) + ^t(P)- 
The term Ec(-P) is the contribution of the 1/q 2 term in the longitudinal piece of the gluon 
propagator. This corresponds to the instantaneous Coulomb interaction term. Explicitly, it 
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takes the form: 

E C (P E ) = g 2 C 2 (F) j {dQ E }^ lo Sf\K E ) lo (C.12) 
Once again, using the spectral representation of the fermion propagator, we obtain: 

/I f°° dk 
(rfq W ^ {i ~ nf{qQ ~ ijL))pa{K)iai0 (c - i3) 

This has no term proportional to g 2 T 2 . Thus, we will not pursue its computation any further, 
but we will rather focus on the transverse part. As for the longitudinal part, p E = at tree 
level, so that it gives a vanishing contribution - This is not true for dressed propagators, and 
in particular, p E does contribute to the quark damping rate. After employing the spectral 
representation of the fermionic propagators, we obtain for the transverse self-energy: 

Zt(Pe) = -9 2 C 2 (F) |(dq) ^(^■-g i g i )( 7 W)po(^K(g ) q)x 

l + n b (q )-n f (k°-fi) 

k° + q° - (iu r + /i) 1 ' 1 

Retarded boundary conditions are obtained from the above by taking the analytic con- 
tinuation iu r +//!—>• p° + ie. Simplifying the Dirac structure and applying the hard thermal 
loop approximation as in the gluon self-energy calculation, we get: 

Er( P) = fa + a / (c.15) 

v ; 8 V ) J An p° - v -p + ie v ; 

with = (1, v). 
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